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SUMMARY 


Laboratory  of  the  National  Aeronautical  Establishment  (NAE)  has  operated 
the  Twin  Otter  atmospheric  research  aircraft  fitted  with  a  CO £1  analyzer 
in  order  to  assess  the  feasibility  of  making  airborne  carbon  dioxide  flux 
measurements.  This  report  will  summarize  progress  to  date  and  present  an 
analyse  of  the  data  collected  during  a  series  of  nine  test  flights  in  August 
and  September,  1982. 

The  experimental  results  show  that  it  is  possible  to  make  realistic 
measurements  of  the  heat  and  C02  fluxes  from  an  aircraft  instrumented  to 
measure  vertical  gusts.  Special  tests  and  analysis  procedures  were  conducted 
to  investigate  the'  accuracy  of  the  carbon  dioxide  measurements.  Recom¬ 
mendations  for  improvements  in  the  design  of  future  test  programs  are 
made. 


RESUME 


% 

A  la  demande  d’Agriculture  Canada,  le  Laboratoire  de  recherche 
en  vol  de  l'Etablissement  aeronautique  national  (EAN)  a  6valu6  la  faisabilite 
de  mesurer  en  vol  le  flux  de  dioxyde  de  carbone  a  l’aide  de  l’avion  de 
recherche  atmospherique  Twin  Otter,  muni  d’un  doseur  de  C02.  Le  rapport 
resume  les  progres  realises  jusqu’ici  et  present*  une  analyse  des  donr.ccs 
obtenues  au  cours  de  neuf  vols  d’esaais  qui  ont  lieu  en  aoQt  et  en  septembre 
1982. 


Les  r&ultats  experimentaux  indiquent  qu’il  est  possible  de  faire 
des  mesures  r&alistes  des  flux  de  chaleur  et  de  C02  a  partir  d’un  avion  muni 
d’instruments  capables  de  mesurer  les  rafales  verticales.  On  a  utilise  des 
m6thodes  speciales  d’essai  et  d ’analyse  pour  ^valuer  la  precision  des  quantity 
de  dioxyde  de  carbone  mesurees.  On  fait  des  recommandations  sur  les 
ameliorations  &  apporter  lore  de  1 ’elaboration  des  programmes  d'essai  future. 
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i.  INTRODUCTION 


The  increasing  level  of  carbon  dioxide  (CO2)  in  the  atmosphere 
has  led  to  much  scientific  speculation  on  long  term  effects  as  well  as 
the  formulation  of  analytical  models  to  predict  the  consequences  on 
global  temperature,  the  weather,  and  the  polar  ice  caps.  Although 
reasonably  accurate  estimates  can  be  made  for  the  input  of  carbon 
dioxide  to  the  atmosphere  from  the  combustion  of  carbon-based  fuels, 
much  less  well  understood  are  the  rates  of  exchange  between  the 
atmosphere  and  the  oceans  and  vegetation. 

For  several  years  the  Agrometeorology  Section  of  Agriculture 
Canada  has  been  developing  instrumentation  and  conducting  ground- 
based  experiments  to  measure  the  CO2  exchange  rate  between  the 
atmosphere  and  various  crops.  In  1979  the  Flight  Research  Laboratory 
of  the  National  Aeronautical  Establishment  (NAE)  was  asked  by 
Agriculture  Canada  (Ref.  1)  to  assess  the  feasibility  of  combining  a 
CO2  sensor  with  the  gust  measuring  capability  of  an  atmospheric 
research  aircraft  to  make  airborne  measurements  of  the  flux  of  carbon 
dioxide  above  various  ecosystems.  If  successful,  an  aircraft  offers  the 
potential  to  quickly  determine  the  CO2  exchange  rate  over  large 
areas  for  inputs  to  theoretical  models,  as  well  as  to  open  the 
possibility  of  inferring  the  growth  rates  of  crops  and  forests  from 
these  flux  measurements. 

In  19S0  the  NAE  Twin  Otter  was  used  in  a  series  of  nine  test 
flights,  five  with  a  CO2  analyzer  built  for  the  Bedford  Institute  of 
Oceanography,  and  four  with  an  open  path  analyzer  developed  by 
Agriculture  Canada.  Although  the  flights  were  restricted  to  the  end 
of  the  growing  season  because  of  other  aircraft  commitments, 
valuable  lessons  were  learned  about  the  use  of  these  sensors  in 
airborne  applications  and  suggestions  were  made  for  improvements 
(Reference  2  and  3).  Resulting  fluxes  calculated  by  the  eddy 
correlation  technique  were  encouraging  and  further  tests  were 
planned. 


In  1981  another  series  of  13  flights  totalling  10  flying  hours  was 
flown  using  only  the  Agriculture  Canada  developed  sensor.  Several 
instrumentation  problems  were  experienced  which  led  to 
recommendations  for  additional  improvements,  but  only  limited 
further  analysis  of  the  flux  data  was  attempted. 

After  undergoing  aerodynamic  and  electronic  refinements,  the 
open  path  CO2  analyzer  was  again  installed  in  the  Twin  Otter  for  a 
ten  hour  series  of  test  flights  in  August  and  September,  1982.  The 
signal  to  noise  ratio  of  the  CO2  analyzer  was  much  improved. 
Consequently,  in  order  to  continue  the  assessment  of  the  use  of 
aircraft  to  make  CO2  flux  measurements,  the  author  has  undertaken 
the  analysis  of  77  of  the  runs  flown  in  this  series  of  tests.  The  results 
are  documented  in  this  report  and  Reference  4,  and  foster  suggestions 
for  further  instrumentation  teste  and  an  improved  experimental  design 
for  the  next  test  series.  A  major  part  of  the  Laboratory's  continued 
interest  in  these  teste  stems  from  the  fact  that,  in  this  particular  test 
environment  (very  low  altitude,  low  winds,  and  small  vertical  gusts 
damped  by  proximity  to  the  surface),  valuable  experience  is  being 
gained  on  probing  the  limits  of  the  capability  of  an  aircraft  to  make 
accurate  gust  and  flux  measurements. 

2.  INSTRUMENTATION 

The  Twin  Otter  measures  the  three  orthogonal  components  of 
atmospheric  motion  in  aircraft-  and  earth-fixed  axes  utilizing  a  nose- 
mounted  gust  boom,  a  3-axis  Doppler  radar,  and  a  strapped-down 
inertial  system  consisting  of  accelerometers,  and  rate  and  attitude 
gyros.  Atmospheric  state  parameters  such  as  temperature  and  dew 
point  are  also  measured,  as  are  the  aircraft  motion,  altitude  (pressure 
and  radio  altimeter),  and  geographical  position  (GNS-500A).  Three  on¬ 
board  microprocessors  compute  these  parameters  in  real  time  and 
present  one-second  averages  to  the  aircrew  on  a  cockpit-mounted 
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plasma  display  unit.  Data  are  processed  by  anti-aliasing  low-pass 
filters  with  a  5  Hz  cutoff,  then  digitized  and  recorded  at  a  rate  of  16 
samples  per  second  in  serial-digital  format  on  two  tracks  of  a  light¬ 
weight  and  reliable  studio  recorder.  Table  1  lists  the  parameters 
recorded  during  these  tests,  track- 1  for  the  raw  sensor  data  and  track- 
2  for  computed  time  histories  output  from  the  microprocessors.  A 
complete  description  of  the  aircraft,  instrumentation,  software,  and 
playback  data  available  is  presented  in  Reference  5. 

The  Agriculture  Canada  open-path  CO 2  analyzer,  based  on  the 
differential  absorption  by  CO2  of  infrared  radiation  at  wavelengths  of 
4.3  and  4.7  pm,  was  mounted  through  the  escape  hatch  in  the  cabin 
roof  of  the  Twin  Otter  (Fig.  1).  With  the  radiation  source  and  detector 
just  inside  the  cabin  and  the  infrared  beam  reflected  by  a  mirror  0.75 
m  above  the  fuselage,  an  in-air  optical  path  length  of  1.5  m  was 
achieved.  An  output  voltage  proportional  to  carbon  dioxide 
concentration  was  fed  to  the  NAE  signal  conditioner  where  it  was 
amplified,  digitized  and  recorded  on  both  tracks  of  the  serial-digital 
recorder.  Since  a  bias  adjustment  was  made  on  the  sensor  after 
attaining  flight  speed  and  the  test  altitude,  the  recorded  signal 
represented  the  fluctuations  in  the  concentration  of  CO2,  with  no 
attempt  made  to  document  the  mean  atmospheric  concentration. 
After  amplification,  the  calibration  of  the  CO2  output  signal  as 
recorded  represented  changes  of  17.54  ppm  per  volt,  giving  a  resolution 
of  0.343  ppm  per  bit  and  a  full  scale  range  of  ±175  ppm. 

The  1980  and  1981  tests  showed  that  the  output  voltage  from  the 
CO2  analyzer  contained  an  unacceptable  level  of  noise,  the  amplitude 
of  which  was  proportional  to  airspeed  (Ref.  2).  This  suggested  that  the 
noise  was  a  result  of  the  inconsistent  reception  by  the  detector  of  the 
reflected  infrared  beam,  due  to  turbulence-generated  vibration  of  the 
mirror  and  its  support  structure.  Prior  to  the  1982  tests,  the  vertical 
tubes  supporting  the  mirror  and  baffles  were  streamlined  (Fig.  2),  and 
noise  levels  appear  to  be  significantly  reduced. 


3 


Incident  solar  radiation  was  measured  by  a  light  sensor  mounted 
vertically  in  the  alt  cabin  roof.  Its  uncaiibrated  signal  was  recorded 
on  both  tracks  of  the  NAE  recorder  and  used  primarily  to  document 
the  passage  of  cloud  shadows.  Another  device  provided  by  Agriculture 
Canada  to  measure  both  incoming  and  reflected  radiation  was  mounted 
on  the  cabin  floor  with  its  sensor  looking  vertically  downward  through 
an  opening  in  the  cabin  floor  and  fuselage.  It  was  used  to  obtain  a 
relative  measurement  of  the  leaf  area  index  of  green  vegetation.  As 
its  signal  was  internally  recorded,  it  was  not  processed  by  NAE. 

3.0  FLIGHTS 

The  aircraft  was  committed  to  another  program  through  3une 
and  July,  so  flights  had  to  take  place  late  in  the  growing  season 
between  August  13  and  September  10.  Nine  project  flights,  totalling 
just  under  ten  flying  hours,  were  conducted.  Table  2  presents  a  brief 
summary  of  the  dates  and  flights  times,  weather  conditions,  and  the 
types  of  terrain  and  vegetation  studied. 

Data  from  the  first  three  flights  were  unuseable  because  of 
instrumentation  problems.  On  Flight  4  several  runs  were  made  over 
Corn  Field  C  near  Russell  Airport,  but  near  the  end  of  the  flight  the 
CO2  analyzer  failed  during  the  only  investigation  of  soybean  fields  in 
the  project.  On  August  26  the  aircraft  was  operated  from  Kingston  to 
make  long  runs  over  Lake  Ontario  to  study  the  noise  level  and 
resolution  of  the  CO2  signal  in  an  area  of  expected  minimal 
fluctuations  in  carbon  dioxide.  This  day  was  almost  overcast  and  had 
the  strongest  winds  of  the  experiment.  On  the  return  flight  to  Ottawa 
several  long  runs  were  made  over  mixed  forest.  On  the  first  part  of 
Flight  7  on  August  31,  two  of  the  Agriculture  Canada  special  test  plots 


near  the  Cedarview  Road  southwest  of  Ottawa  were  overflown.  The 
first  of  these^  Com  Field  A,  was  the  same  field  over  which  ground- 
based  flux  measurements  were  made  prior  to  the  flight  tests,  utilizing 
the  seme  CO2  analyzer  carried  aboard  the  aircraft.  The  second,  Corn 
Field  B,  was  centred  about  one  kilometer  northwest  of  A.  Flight  7 
concluded  with  a  pair  of  long  reciprocal  passes  over  the  Larose  Forest 
and  Mer  Bieue,  a  marshy  area,  southeast  of  Ottawa. 

The  final  two  flights,  although  in  September,  had  the  clearest 
skies  of  the  series  and  temperatures  typical  of  the  earlier  flights.  The 
winds  were  very  light,  especially  on  Flight  8.  On  September  8,  five 
fairly  long  runs  were  performed  over  mixed  farmland  near  Richmond, 
Ontario,  followed  by  several  passes  over  Corn  Field  A  and  Lac 
Deschenes.  On  September  10,  eight  runs  each  were  made  over  Corn 
Fields  A  and  B  and  a  mixed  forest  just  west  of  Field  B. 

Since  the  1980  and  1981  studies  showed  that  the  CO2  sensor  had 
noise  levels  directly  related  to  airspeed,  the  1982  tests  were  flown  at 
the  lowest  possible  safe  airspeed,  utilizing  about  15  degrees  of  flap. 
On  most  of  the  com  field,  farmland,  and  forest  runs,  true  airspeeds 
averaged  about  50  m  s*l  .  For  safety  reasons  the  Lake  Ontario  series 
was  conducted  at  higher  airspeeds,  around  60  m  s~l  .  Care  was  taken 
to  minimize  altitude  excursions.  Very  low  altitudes  were  achieved  on 
some  runs,  averaging  less  than  8  m  on  the  September  10  com  field 
runs,  for  example.  Greater  altitude  variations  were  unavoidable  over 
the  Richmond  farmland  and  the  forests  because  of  the  presence  of 
power  lines  and  tall  trees. 

*.0  DATA  ANALYSIS 

Prior  to  the  summer  series  of  project  flights,  several 
improvements  were  made  to  the  Twin  Otter  playback  software.  On  a 
routine  basis  after  each  flight,  the  real-time  computed  data  recorded 


on  track-2  are  used  to  produce  a  flight  path  plot  and  a  printout  of 
winds,  temperatures,  positions  etc.  For  the  event-marker-in  portions, 
which  denote  particular  measurement  runs,  RMS  and  peak  gust 
velocities  are  computed  along  with  run  averages  of  all  pertinent 
parameters.  For  these  tests  the  program  was  modified  to  include  the 
run-average  CO2  vertical  flux  computed  from  the  relationship: 


pwco2 
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Where  Wge  is  the  vertical  component  of  air  velocity,  CO2  is  the 
carbon  dioxide  concentration,  and  P  is  the  density  of  air.  The  overbars 
represent  run  averages. 

More  detailed  analysis  of  the  flight  data  takes  place  on  the 
NRC  Time-Shared  IBM  3033  System,  where  many  programs  and 
subroutines  have  been  developed  to  utilize  data  recorded  on  track- 1  of 
the  Twin  Otter  recorder  for  the  computations  of  gusts,  fluxes,  power 
spectra,  etc.  New  software  has  now  been  added  that  uses  the 
identical  subroutines  to  produce  the  same  outputs,  but  which  employs 
the  track-2  airplane-computed  data  as  inputs. 

For  the  1982  test  data,  it  was  initially  planned  to  perform  all 
the  analysis  using  only  the  track-2  parameters.  It  was  discovered, 
however,  that  on  many  of  the  runs,  the  real-time  computed  wind 
components  possessed  biases  at  the  start  that  took  several  seconds  to 
decay.  The  resultant  drifts  are  a  product  of  the  low  filter  breakpoint 
frequency  (0.02HZ)  required  in  the  real-time  complementary  digital 
filtering  routine  used  to  compute  the  aircraft  inertial  velocity 
components.  The  direct  cause  is  the  tight  manoeuvring  needed  to 
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position  the  aircraft  at  low  altitude  over  the  usually  small  fields  while 
avoiding  obstacles  at  their  margins.  Any  false  drift  in  the  vertical 
gust  component  clearly  leads  to  inaccurate  carbon  dioxide  and  heat 
fluxes,  especially  for  the  short  corn  field  runs  where  the  drift  often 
existed  for  up  to  a  third  of  the  recorded  run.  Consequently,  it  was 
decided  to  analyze  all  the  runs  using  the  raw  track- 1  data  and  to 
recompute  the  gusts  and  fluxes  on  the  NRC  Time-Shared  System  (TSS). 
Although  the  complementary  filtering  technique  used  in  these 
programs  is  virtually  identical  to  that  in  the  aircraft  real-time 
processing,  drifts  can  be  avoided  by  initiating  the  TSS  analysis  at  the 
point  where  the  aircraft  has  been  stabilized  over  the  test  field. 

The  three  components  of  the  true  gust  velocity  in  airplane-  and 
earth-fixed  axes  are  computed  by  the  methods  detailed  in  Reference 
5.  For  these  data  the  earth  axis  system  used  has  the  Uge  axis  aligned 
north/south  with  a  positive  gust  from  the  north.  A  positive  lateral 
component,  Vge,  is  from  the  east,  and  the  positive  sense  for  the 
vertical  gust  is  upwards.  Prior  to  the  CO2  flux  calculation,  the  run 
average  was  removed  from  the  vertical  gust  and  CO2  signals.  The 
mounting  location  for  the  CO2  analyzer  was  6.5  m  aft  of  the  angle  of 
attack  sensor  on  the  gust  boom.  In  the  CO2  flux  calculation, 
therefore,  the  Wge  and  CO2  signals  were  time  correlated  by 
advancing  the  CO2  signal  by  the  integral  number  of  data  slices 
(usually  2)  closest  to  6.5*16/(true  airspeed). 

Run-average  heat  and  CO2  fluxes  were  calculated  by  two 
procedures.  Method  A  was  to  simply  average  the  product  over  the 
entire  event-in  portion  of  the  run,  i.e. 

_  N  (Wge  (i)-Wge)  *  (CO  (i)-COT)  *  p 

pvco2  .  f=1  - i - 
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In  Method  B,  the  co-spectrui  *>f  the  product  was  integrated  over  a 
selected  wavelength  range  from  10  to  800  m.  Since  the  project  runs 
varied  in  length  from  about  one  to  greater  than  13  km,  a  meaningful 
comparison  of  fluxes  between  runs  over  different  terrains  can  really 
only  be  achieved  by  using  a  common  limited  wavelength  band.  It  is 
important  to  note  that  all  of  the  RMS  values  and  heat  fluxes  presented 
in  this  report  were  computed  by  Method  B.  To  show  the  above- 
mentioned  variations,  three  computations  of  the  carbon  dioxide  flux 
are  given  for  each  run,  i.e.  track-2  and  Methods  A  and  B  for  track-1 
(Tables  3-7).  For  the  final  summaries  (Tables  8-10),  only  fluxes 
calculated  over  the  10  to  800  m  wavelength  range  are  included.  No 
attempt  was  made  to  correct  the  CO2  fluxes  for  the  effects  of  the 
sensible  and  latent  heat  fluxes,  so  the  resultant  CO2  fluxes  could  be 
overestimated  by  up  to  twenty  per  cent  over  land  and  possibly  more 
over  water  (Ref.  6). 

A  total  of  77  runs  was  analyzed.  Outputs  included  printed 
listings  of  winds,  temperatures,  flight  conditions,  RMS  gust  velocities, 
flux  estimates,  spectra,  and  co-spectra.  Time  histories  of  the 
following  10  parameters  were  written  to  disk  and  backed-up  to 
magnetic  tape  for  storage:  event  marker,  3-axis  gust  components, 
radiometer,  radio  altimeter  height,  temperature,  CO2  signal,  CO2 
flux,  and  the  cumulative  average  of  the  CO2  flux.  The  last  parameter 
was  computed  from  the  relationship: 

SUMF  (i)  *  E_1  Flux(j)  Where  Flux(j )=(Wge( J )  -  Wge)«(C02(j )  -  CC>2)«p 


Time  history  plots  of  the  important  parameters  were  prepared 
for  all  the  runs,  examples  of  which  are  shown  in  Figures  3  to  8  for  the 
various  terrains  studied.  It  must  be  noted  that  the  radio  altimeter 
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height  shown  (RALT)  is  measured  from  the  bottom  of  the  fuselage  to 
the  surface.  Since  the  CO2  sensor  is  mounted  on  the  top  of  the 
fuselage,  the  actual  carbon  dioxide  sample  was  taken  2.3  m  above  the 
height  shown  in  the  figures  and  tables  below.  Data  from  only  70  of  the 
77  runs  are  presented  in  this  report.  There  were  occasional  false 
spikes  in  the  vertical  gust  and  CO2  signals,  such  as  seen  in  Figure  8. 
Examination  of  the  flux  and  SUMF  time  histories  showed  that  on  most 
runs  these  made  no  significant  contribution  to  the  flux,  but  where  the 
opposite  was  true,  the  run  was  excluded  from  further  analysis. 

Gust,  CO2  ,  and  temperature  spectra,  and  the  Wge*CC>2  and 
heat  co-spectra  were  output  to  a  labelled  dataset  for  each  run.  A 
special  program  was  then  used  to  produce  average  unnormalized 
spectra  and  co-spectra  for  a  particular  series  of  runs,  e.g.  for  a 
particular  terrain  type  (Figures  9  to  18). 

Copies  of  the  time  histories  on  magnetic  tape,  analog  plots,  and 
the  summary  tables  have  been  forwarded  to  Macdonald  College  of 
McGill  University  for  their  analysis  under  contract  to  Agriculture 
Canada. 

5.0  DISCUSSION  OF  RESULTS 


Tables  3  to  7  present  the  computed  results  for  each  run, 
grouped  by  terrain  type  and  flight.  After  each  series  of  runs  on  a 
given  flight,  the  average  values  for  all  runs  are  given  along  with  the 
standard  deviation,  o  ,  to  show  the  variation  between  runs.  These 
averages  and  standard  deviations  are  carried  to  the  summary  Table  8 
to  facilitate  comparison  of  results. 
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Perhaps  the  most  noticeable  observation  from  Tables  3  to  7  is 
the  large  run-to-run  variation  in  the  measured  CO2  flux  for  all  three 
methods  of  computation.  In  most  cases  the  track-2  CO2  fluxes 
display  greater  variability  than  the  track- 1  results,  principally  due  to 
the  drifts  in  the  track-2  vertical  wind  component  on  some  runs  as 
discussed  in  Section  4.0.  As  would  be  expected,  the  average  track- IB 
fluxes  are  generally  smaller  than  those  computed  from  track- 1  data 
using  Method  A,  since  the  contributions  by  eddies  with  wavelengths 
longer  than  800  m  are  excluded.  The  heat  flux  data  appear  to  be  less 
variable  than  the  CO2  fluxes.  The  RMS  of  the  vertical  gust  velocity 
and  the  CO2  concentration  tend  to  be  fairly  consistent  from  run  to 
run.  Throughout  all  the  tests,  the  standard  deviation  of  the  RMS 
vertical  gust  velocity  averages  only  about  0.1  m  s"f,  suggesting  that 
the  considerable  variability  in  the  fluxes  is  not  related  to  the  vertical 
gust  component. 

Despite  the  variability  of  the  fluxes  between  runs,  the  results 
of  this  feasibility  study  should  be  considered  encouraging.  The 
average  heat  and  carbon  dioxide  fluxes  appearing  in  Table  8  are  of  the 
correct  sign  and  of  reasonable  magnitude,  comparing  favourably  with 
previous  ground-based  measurements.  In  the  latter  case  it  has  been  an 
accepted  practice  to  average  the  results  of  several  runs,  as  has  been 
done  here.  This  appears  to  be  especially  necessary  when  only  short 
runs  can  be  made,  as  was  the  case  for  the  corn  fields. 

The  analysis  cf  flux  data  has  usually  been  based  on  the 
assumption  of  steady  state,  i.e.  a  time  invariant  boundary  layer  in 
equilibrium  with  the  sources  and  sinks  of  CO2  and  heat.  This  is 
unlikely  to  have  been  the  case  for  these  tests,  for  wind  shifts,  moving 
cloud  shadows,  etc.  can  cause  singular  events  leading  to  the  over  or 
underestimation  of  the  turbulent  fluxes.  For  small  fields  with 
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marginal  fetch,  wind  shifts  can  bring  air  from  adjacent  fields  with 
different  vegetation  and  presumably  different  CO2  concentrations. 
Although  there  are  only  a  limited  number  of  long  runs  in  the  1982  study 
(Larose,  Richmond,  and  Lake  Ontario),  these  data  and  those  of 
Reference  3  support  the  observation  of  improved  consistency  with  the 
length  of  the  flight  path.  Examination  of  the  traces  of  the  cumulative 
average  CO2  flux  (SUMF)  for  the  longer  runs  reveals  that  at  least  30 
seconds,  and  often  a  minute  (3  km),  is  required  for  them  to  reach  a 
stable  level  representative  of  their  final  value. 

Run  04-08C  is  a  good  example  of  a  singular  event  discussed 
above,  one  with  the  largest  track-lB  CO2  flux  and  second  largest  heat 
flux  of  the  entire  experiment  (Table  4,  Fig.  4).  In  the  26-30  second 
range  of  the  time  histories,  there  is  a  200  m  wide  updraft  averaging 
about  2  m  s-l  ,  well  correlated  with  an  increase  in  temperature  and  a 
decrease  in  CO2  concentration.  None  of  the  other  runs  over  Corn 
Field  C  that  day  showed  such  an  event,  nor  can  anything  unique  be 
found  in  the  conditions  for  this  particular  run.  It  should  be  noted, 
however,  that  the  average  CO2  flux  for  all  the  passes  over  this  field 
(-46  kg  ha-1  h_1  ,  Table  8)  is  the  largest  recorded  in  these  tests.  This 
is  the  earliest  flight  of  those  reported  here  (August  23)  and  would  be 
expected  to  exhibit  the  largest  CO2  exchange  rates.  Also  the  winds 
averaged  209  degrees  at  6.7  m  s-l  ,  the  strongest  for  all  the  corn  field 
runs,  and  probably  contributed  to  the  increased  flux  levels. 

Figures  9  to  17  depict  average  co-spectra  for  the  heat  and  CO2 
fluxes  for  the  various  types  of  terrain  studied.  The  standard  deviation 
of  the  spectral  estimates  at  each  frequency  is  relatively  large  since 
the  variability  of  the  fluxes  themselves  is  considerable.  Nevertheless, 
these  presentations  are  useful  for  showing  the  size  of  the  eddies 
contributing  to  the  fluxes  and  for  examining  the  potential  differences 
between  vegetation  types. 
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At  wavelengths  shorter  than  about  30  m,  there  is  a  negligible 
contribution  to  the  CO2  flux,  but  a  small  but  significant  contribution  to 
the  heat  flux.  This  is  an  indication  that,  for  an  airspeed  of  50  m  s~l,  the 
upper  limit  of  the  frequency  response  of  the  CO2  analyzer  is  about  1  to 
1-1/2  Hz  in  the  airborne  application.  This  is  lower  than  its  reported 
performance  in  ground-based  studies,  and  suggests  the  need  for 
additional  tests  of  this  aspect  of  its  operation. 

For  most  of  the  plots  the  heat  and  CO2  co-spectra  show  a 
similar  dependence  on  wave  number,  but  are  opposite  in  sign.  This 
indicates  that  basically  the  same  eddies  are  responsible  for  the 
vertical  transfer  of  heat  and  carbon  dioxide,  and  that  during  the 
photosynthetically  active  part  of  the  day,  the  vegetation  is  a  source  of 
heat  and  a  sink  for  CO2  • 

The  co-spectra  show  that  the  bulk  of  the  CO2  and  heat 
transfer  is  carried  by  eddies  with  wavelengths  in  the  50  to  500  m 
range.  There  are  observable  differences  in  the  locations  of  the  peaks 
for  different  days  and  terrains,  which  tend  to  be  more  evident  in  the 
heat  co-spectra  because  they  are  less  variable.  A  comparison  of 
Figures  9  and  11,  for  example,  shows  heat  and  CO2  transfer  at  longer 
wavelengths  over  Field  C  than  Field  A.  If  this  trend  is  real,  it  could 
be  due  to  the  higher  altitude  flown  over  C,  the  stronger  winds,  or 
perhaps  the  greater  photosynthetic  activity  earlier  in  the  season. 
Similarly  the  runs  over  the  Larose  Forest  illustrate  an  increased 
contribution  to  the  heat  flux  by  long  wave  length  eddies  up  to  800  m 
(Fig.  13).  Examination  of  the  vertical  temperature  profiles  recorded 
by  the  aircraft  to  a  height  of  300  m  shows  the  greatest  gradients  on 
the  day  Com  Field  C  was  studied,  while  the  Larose  runs  were  flown  on 
the  day  with  the  second  steepest  lapse  rate.  Reference  7  indicates 
that  a  shift  to  longer  wavelengths  in  the  co-spectra  can  be  expected 
with  increased  atmospheric  instability. 
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Confidence  in  the  results  and  in  the  use  of  the  eddy  correlation 
technique  to  compute  the  fluxes  is  bolstered  somewhat  by  the  co¬ 
spectra  averaged  for  the  10  runs  over  water  (Fig.  12).  As  would  be 
anticipated,  the  co-spectral  estimates  and  mean  fluxes  over  water  are 
an  order  of  magnitude  lower  than  those  for  the  runs  over  vegetation. 
There  is  no  noticeable  preference  for  contributions  at  any  particular 
wavelength,  indicating  an  absence  of  correlated  error  between  the 
vertical  gust  component  and  the  C02  or  temperature  signals. 

Tables  9  and  10  present  the  final  summary  of  the  average 
carbon  dioxide  and  heat  fluxes  versus  flight  number  and  terrain  type. 
It  is  difficult  to  draw  firm  conclusions  from  these  data.  There  are  an 
insufficient  number  of  entries  in  each  category  to  positively  relate 
observed  variations  to  physical  explanations.  Trends  suggested  by  the 
data  in  one  category  are  sometimes  contradicted  by  the  limited  results 
from  another  category.  Obviously  there  are  many  factors  that  can 
affect  the  heat  and  carbon  dioxide  fluxes.  These  include  the  altitude 
flown,  the  incident  radiation,  wind  speed  and  direction,  the 
temperature  and  moisture  content  of  the  soil,  plant  stress, 
atmospheric  stability,  and  the  phase  of  the  growing  season.  As  this 
program  was  basically  an  instrument  test  and  feasibility  study,  not  all 
of  these  factors  were  documented  to  the  degree  required  for 
meaningful  correlation  with  the  average  fluxes.  It  is  planned  to  do  so 
in  the  next  test  program.  It  is  also  recommended  that,  in  the  next 
CO2  flux  measuring  program,  the  test  areas  be  restricted  to  a  number 
that  can  be  investigated  on  every  flight.  With  entries  at  all  possible 
array  locations  in  summary  tables  such  as  9  and  10,  significant  trends 
in  the  data  may  then  become  apparent. 


6.0  SPECIAL  TESTS 


An  important  aspect  of  this  feasibility  study  was  the 
investigation  of  possible  contributors  to  the  run  to  run  flux  variations. 
At  this  point  the  subject  of  the  reliability  and  accuracy  of  the  CO2 
analyzer  should  be  addressed  in  more  detail.  Particular  questions  to 
be  considered  are  whether  the  streamlining  of  the  mirror  supports  has 
completely  eliminated  noise  problems,  and  whether  the  excursions  in 
the  CO2  analyzer  output  represent  actual  CO2  fluctuations  or  are 
partially  a  result  of  some  other  interference  effect. 

Special  tests  were  performed  to  attempt  to  answer  these 
questions.  These  included  low  altitude  runs  over  water  distant  from 
shore  where  real  CO2  fluctuations  would  be  expected  to  be  a 
minimum.  The  results  in  Table  5  are  promising,  for  they  show  the 
anticipated  very  low  heat  and  CO2  fluxes  with  small  standard 
deviations,  especially  when  compared  with  the  over-land  runs  (Table 
8).  However,  the  eddy  correlation  technique  used  to  compute  the 
fluxes  can  mask  noise  problems.  This  is  one  of  its  advantages  under 
normal  circumstances.  A  signal  contaminated  with  random  noise  can 
be  multiplied  by  another  to  give  a  good  average  flux  estimate,  as  long 
as  the  noise  on  the  one  signal  is  uncorrelated  with  the  second  signal. 
The  flight  results  therefore  indicate  that  the  measured  CO2 
concentration  fluctuations  are  either  real  or,  if  they  have  a  false 
component,  it  is  uncorrelated  with  the  vertical  gust  velocity,  so 
realistic  flux  estimates  still  result. 

6.1  Frequency  Response  of  the  CO2  Analyzer 

Examination  of  Figures  3  and  4  reveals  that  the  CO2  trace  has 
less  high  frequency  content  than  the  gust  and  temperature  time 
histories.  To  further  investigate  this  point,  average  power  spectra  for 
the  vertical  gust,  temperature,  and  CO2  signals  were  computed  for  23 
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runs  over  Corn  Field  A  (Figure  18).  The  Wge  and  temperature  spectra 
show  the  expected  -5/3  slope  at  high  frequencies;  the  CO2  spectra  for 
the  corn  field  runs  and  all  forest  runs  do  not.  The  steep  slope  at  k  >  0.02m- 1 
indicates  a  reduced  response  at  frequencies  greater  than  about  1  Hz. 

If  turbulence  generated  by  the  flow  about  the  mirror  supports  and 
baffles  caused  mirror  vibration  and  resultant  noise  in  the  CO2  output 
signal,  it  would  be  expected  to  appear  at  these  frequencies.  The  lack 
of  high  frequency  fluctuations  could  be  an  indication  that  the 
streamlining  of  the  mirror  supports  has  had  the  desired  effect  of 
improving  the  signal  to  noise  ratio.  The  steep  slope  of  the  CO2 
spectrum  is  of  concern,  however,  and  deserves  additional  investigation 
in  future  tests. 

6.2  CO2  Fluctuations  and  Sideslip  Angle 

The  average  RMS  value  for  the  CO2  concentration  is  given  in 
Table  8  for  each  series  of  runs.  Over  land  the  values  are  from  3  to  4 
ppm.  Over  water  they  are  2.3  to  2.4  ppm,  higher  than  expected  and 
perhaps  an  indicator  of  a  level  of  error.  Close  examination  of  the  CO2 
traces  shows  that  many  of  the  signal  fluctuations  occur  at  frequencies 
near  the  short  period  and  lateral  oscillatory  modes  of  the  aircraft  (0.3 
to  0.5  Hz).  Streamlining  the  mirror  supports  might  have  made  the 
analyzer  sensitive  to  lateral  forces  caused  by  sideslip.  If  this  were  the 
case,  the  CO2  signal  would  be  expected  to  increase  with  either 
positive  or  negative  sideslip. 

The  analysis  software  was  modified  to  compute  sideslip  angle  p 
at  the  mounting  location  of  the  CO2  analyzer.  The  correlations 
between  the  CO2  signal  and  p  the  lateral  gust  velocity  were  also 
computed.  Three  runs  were  analyzed  with  this  modified  software,  Run05-11L 
over  Lake  Ontario,  Run  07-01F  over  the  Larose  Forest,  and  Run  09-CI, 
a  special  test  in  which  the  aircraft  was  flown  in  relatively 
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smooth  air  at  1200  ft  above  ground  and  subjected  to  pitch,  roll,  and 
yaw  control  inputs.  Figure  19  shows  the  time  histories  for  the  three 
runs.  The  top  two  traces  represent  the  lateral  and  vertical  gust 
velocities  computed  in  earth  axes  with  the  x  axis  aligned  with  the  run- 
averaged  heading  of  the  aircraft.  The  Vge  trace  approximates  the 
lateral  gust  velocity  experienced  by  the  CO2  sensor.  The  third  and 
fourth  traces  depict  the  sideslip  angle  (3  at  the  sensor  and  the  CO2 
signal  respectively.  It  can  be  seen  that  substantial  sideslip  angles  are 
generated,  up  to  about  7  deg.  at  'A'  in  the  forest  run  and  after  'B' 
during  the  yawing  manoeuvres  in  the  control  input  tests.  Although  at 
some  points  there  is  the  suggestion  of  a  relationship  between  the 
sideslip  angle  and  the  CO2  signal,  there  does  not  appear  to  be  a 
conclusive  one  to  one  correspondence  between  p  excursions  and 
positive  CO2  fluctuations.  Also,  in  the  control  input  case,  there  is  no 
noticeable  increase  in  CO2  activity  after  'B'  when  the  yawing 
manoeuvres  commence. 

The  bottom  three  traces  in  Figure  19  show  the  products  of  CO2 
times  p,  lateral  gust  velocity,  and  vertical  gust  velocity  respectively. 
The  lateral  flux  of  CO2  tends  to  show  larger  excursions  than  the 
vertical  flux  for  the  forest  and  lake  runs,  but  about  the  same  for  the 
control  input  case.  This  is  not  unexpected,  for  at  the  low  altitudes 
flown  in  the  first  two  cases,  the  vertical  gusts  are  damped  because  of 
proximity  to  the  surface.  Moreover,  since  each  run  had  a  considerable 
crosswind  component,  the  lateral  flux  of  CO2  should  equal  or  exceed 
the  vertical  transport. 

Since  it  was  postulated  that  CO2  fluctuations  could  be 
correlated  with  sideslip  at  particular  frequencies,  power  spectra  for  p 
and  CO2  have  been  computed  for  the  three  special  test  runs.  The 
spectra  at  the  bottom  of  Figure  20  show  peaks  clearly  associated  with 
the  yawing  oscillations  of  the  aircraft  at  wave  numbers  ranging  from 
3.5  to  6.8  x  10"3  m-l  (frequencies  of  0.19  to  0.32  Hz).  This  is 
especially  evident  for  the  forced  yaws  in  the  control  input  case,  where 


the  peak  in  the  spectrum  is  more  than  an  order  of  magnitude  higher 
than  it  would  have  been  without  the  rudder  inputs.  If  side  forces  cause 
mirror  deflections  that  shift  the  reflected  infrared  beam  partially  off 
the  detector,  then  a  false  positive  CO2  reading  would  result  for  both 
left  and  right  sideslips,  i.e.  the  CO2  spectrum  would  show  a  peak  at 
twice  the  frequency  of  the  p  oscillations.  The  CO2  spectra  in  the 
upper  part  of.  Figure  20  show  no  obvious  peaks  at  frequencies 
corresponding  with  either  the  P  spectral  peak  or  twice  its  wave 
number.  It  can  be  concluded  that  very  little,  if  any,  of  the  CO2  signal 
fluctuation  is  caused  by  deflection  of  the  mirror  support  structure  due 
to  aerodynamic  side  forces. 

6.3  Effects  of  Aircraft  Wake  and  Exhaust  During  Repeated  Runs 

For  the  smaller  fields,  runs  on  reciprocal  headings  have  been 
made  at  intervals  as  short  as  1-1/2  minutes,  bringing  into 
consideration  the  possible  contamination  of  the  data  as  a  result  of  the 
previous  pass  of  the  aircraft.  In  light  wind  conditions,  or  for  runs 
nearly  along  the  wind  direction,  it  is  possible  that  the  atmosphere 
disturbed  by  the  aircraft  wake  will  not  be  returned  to  a  natural 
equilibrium  prior  to  the  next  passage  of  the  aircraft.  Resultant 
disturbances  would  be  difficult  to  detect  and  impossible  to  correct  in 
the  data. 


A  more  serious  potential  problem  with  repeat  passes  is 
contamination  of  the  CO2  readings  by  carbon  dioxide  emitted  in  the 
engine  exhaust  from  previous  runs.  An  example  calculation  will 
illustrate  this.  Assume  that  in  the  1-1/2  minutes  between  passes,  the 
exhaust  released  in  one  second  is  uniformly  mixed  in  a  volume  of  air 
defined  by  3  times  the  aircraft  wingspan  (3  x  20  m),  the  run  altitude 
(12  m),  and  the  distance  flown  in  one  second  (50  m).  At  the  power 
setting  required  for  100-knot  cruise  near  sea  level,  each  PT6  engine 
bums  at  least  230  lb  per  hour  of  fuel  to  produce  450  lb  of  CO2  (Ref. 


8).  In  one  second,  then,  2*  450  *454  =  113  grams  of  CO2  are  produced 

3600 


and  the  resulting  concentration  is: 


113 


g  CO, 


50*12*60 


m' 


cm 


m' 


3  3 

100 3  cnT 


1.2*10“3  g  air 


=  2.6*10 


-6 


2.6*10"6  * 


.-6 


1.7*10”  parts  C0„  per  part  air 


This  1.7  ppm  concentration  is  of  the  same  order  as  the 
measured  RMS  of  the  CO2  signal.  Even  relaxing  the  above 
assumptions  to  allow  a  greater  dissipation  of  the  wake  between  runs 
will  not  eliminate  the  potential  for  contamination  of  the  CO2 
measurements  by  the  exhaust,  for  there  may  be  isolated  pockets  of 
higher  exhaust  concentration  caused  by  some  environmental  condition 
or  the  wingtip  trailing  vortices.  For  the  data  in  this  report,  the 
problem  could  likely  only  occur  during  the  corn  field  runs. 
Fortunately,  in  most  cases  there  was  a  cross-track  wind  component, 
which  even  in  the  worst  case,  would  have  transported  the  aircraft 
wake  240  m  laterally  between  runs.  For  Runs  3A,  5A,  and  7A  of  Flight 
9,  a  diagonal  pass  across  the  field  was  made  just  80  seconds  after  a 
north-south  run.  Although  the  aircraft  wake  was  probably  crossed  on 
these  runs,  there  are  no  apparent  anomoiies  in  the  temperature,  gust, 
or  CO2  concentration  data  that  can  be  conclusively  associated  with 
interception  of  the  wake  or  exhaust.  Nevertheless,  since  the 
possibility  of  data  contamination  does  exist,  it  is  highly  recommended 
that  in  future  CO2  flux  measurements,  longer  fields  should  be  chosen 
and  repeated  passes  avoided.  If  this  is  impracticable,  all  runs  should 
be  flown  crosswind  with  perhaps  at  least  3  minutes  between  passes. 


7.0  CONCLUSIONS 

The  feasibility  of  measuring  carbon  dioxide  fluxes  from  an 
aircraft  has  been  demonstrated.  Carbon  dioxide  and  heat  fluxes 
averaged  for  several  runs  were  of  the  correct  sign  and  had  magnitudes 
consistent  with  previous  ground-based  measurements. 


g  co2 

g  air 


18 


The  aerodynamic  improvements  made  to  the  Agriculture 
Canada  open  path  CO2  analyzer  appear  to  have  improved  its  signal  to 
noise  ratio  by  eliminating  high-frequency,  airspeed-dependent 
fluctuations  associated  with  vibration  of  the  mirror  support  structure. 
Further  flight  testing  of  the  analyzer  is  required  to  determine  why  its 
response  at  frequencies  above  1  Hz  appears  limited  in  the  airborne 
application  compared  with  its  ground-based  performance. 

The  output  signed  from  the  CO2  sensor  possessed  fluctuations 
that  were  suspected  of  being  too  large,  particularly  in  over-water 
tests.  Special  analysis  procedures  demonstrated  that  these  excursions 
were  not  correlated  with  sideslip  angle,  and  were  therefore  not  a 
result  of  lateral  aerodynamic  forces  on  the  mirror  support  structure. 
No  other  correlations  were  found  that  suggested  these  CO2  readings 
were  anything  but  real.  Nevertheless,  it  is  felt  that  a  CO2  analyzer 
should  be  developed  that  is  mounted  entirely  within  the  aircraft 
fuselage  and  samples  air  ducted  to  it  from  an  external  inlet. 

The  computed  heat  and  CO2  fluxes  showed  considerable 
variation  from  run  to  run.  Even  when  averaged  over  all  runs  for  a 
particular  flight  and  vegetation  type,  it  was  difficult  to  relate  the 
resultant  fluxes  to  environmental  conditions,  terrain,  etc.  To  that 
end,  considerable  attention  must  be  given  to  the  design  of  future 
experiments  in  order  to  make  such  conclusions  possible.  The  next  test 
series  should  cover  a  larger  part  of  the  growing  season  to  provide 
greater  variations  in  photosynthetic  activity  (eg,  mid-June  to  the  end 
of  August).  The  number  of  test  areas  (terrains)  should  be  restricted  to 
about  four,  with  each  flown  on  every  flight.  One  of  these  could  be 
adjacent  to  the  NAE  Mobile  Meteorological  Tower,  which  would 
provide  useful  supporting  wind  and  temperature  data  at  four  levels  in 
the  lower  30  m  of  the  atmosphere.  Run  lengths  must  be  at  least  three 
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kilometers,  and  repeated  passes,  if  done  at  all,  must  be  arranged  so  as 
to  avoid  contamination  by  the  turbulent  wake  and  exhaust  of  the 
previous  pass.  For  every  flight  a  sounding  to  about  5000  ft  should  be 
flown  to  document  the  stability  of  the  atmosphere.  Then,  given  a 
sufficient  number  of  test  cases,  average  fluxes  appearing  in  tables 
such  as  9  to  10  in  this  report,  can  be  examined  for  trends  that  can  be 
related  known  physical  processes  and  environmental  conditions. 
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DATA 

RECORDED  ON  CF-POK  -  AGH 

ICULTURAL  PROJECT  -  1982 

Channel 

Track  1 

Track  2 

0 

GMT 

GMT 

1 

9-Level  Event  Marker 

9-Level  Event  Marker 

2 

Doppler  Vx 

GMT  -  hour 

3 

Doppler  Vy 

GMT  -  minute 

4 

Doppler  Vz 

GMT  -  second 

5 

Dynamic  Pressure  (Airspeed) 

Latitude  -  degrees 

6 

Magnetic  Heading 

Latitude  -  minutes 

7 

Radio  Altimeter  Height 

Longitude  -  degrees 

8 

Total  Temperature 

Longitude  -  minutes 

9 

Dew  Point 

Range  from  Reference 

10 

GNS  Latitute 

Bearing  from  Reference 

11 

GNS  Longitude 

True  Heading 

12 

Vertical  Acceleration 

Wind  Speed 

13 

Lateral  Acceleration 

Wind  Direction 

14 

Longitudinal  Acceleration 

Static  Pressure  (corrected) 

15 

Wind  Direction 

True  Airspeed 

16 

Wind  Speed 

Pitch  Attitude 

17 

Total  Temperature  (fine) 

Roll  Attitude 

18 

Angle  of  Attack 

Uge  -  North/South  Wind  Component 

19 

Angle  of  Sideslip 

Vge  -  East /West  Wind  Component 

20 

Pitch  Attitude 

Wge  -  Vertical  Wind  Component 

21 

Roll  Attitude 

*  Vertical  Wind  by  Pressure  Method 

22 

Pitch  Rate 

Static  Temperature 

23 

Yaw  Rate 

Reverse  Flow  Temperature 

24 

Roll  Rate 

Dew  Point 

25 

— 

*  Rate  of  Change  of  Pressure  Height 

26 

C02  Analyzer 

*  Waire 

27 

Radiometer 

Radio  Altimeter  Height 

28 

Static  Pressure 

*  PS25 

29 

True  Airspeed 

C02  Analyzer 

30 

Time  -  Seconds 

Radiometer 

Computed  parameters  recorded  for  special  software  tests  done  on  same  flights. 


At  Ottawa  Airport  at  takeoff. 


TABLE  3 


FLUX  SUMMARY  -  CORN  FIELD  A 


Fit 

Run 

Alt 

m 

Dist 

km 

aWge 

mS1 

ffC02 

ppm 

Heat 

Flux 

Watts  m'2 

co2 

Flux  kg  ha  1  h‘* 

Trk  2 

Trk  1A 

Trk  IB 

n 

mm 

n 

1.63 

0.62 

4.0 

81 

-12 

-11 

-7 

mm 

1.64 

0.53 

3.7 

61 

57 

-11 

-22 

1.80 

0.71 

3.9 

36 

12 

-2 

-6 

4A 

15 

1.66 

0.76 

4.2 

61 

-2 

5 

6 

BS 

5A 

13 

1.01 

0.53 

3.2 

54 

-15 

-8 

-2 

AVG 

16 

1.55 

0.63 

3.8 

59 

8 

-6 

-6 

a 

0.09 

0.3 

14 

29 

6 

9 

8 

■g 

27 

1.54 

0.66 

3.7 

72 

-53 

-35 

-35 

25 

1.62 

0.81 

3.7 

129 

-57 

-58 

-74 

»il 

26 

2.03 

0.63 

3.5 

88 

-5 

-17 

5 

4A 

27 

1.97 

0.76 

3.3 

76 

-32 

-34 

-18 

6A 

25 

1.97 

0.87 

2.8 

93 

WBm 

-17 

-11 

7A 

24 

2.22 

0.61 

3.2 

73 

-14 

-24 

8A 

25 

2.25 

0.87 

2.6 

109 

-19 

-20 

-23 

9A 

27 

2.13 

0.58 

3.5 

65 

-35 

-43 

-38 

10A 

26 

2.16 

0.67 

3.3 

67 

-19 

-19 

11A 

25 

2.17 

0.52 

3.2 

66 

-12 

-12 

mm 

AVG 

26 

2.01 

0.70 

3.3 

84 

-23 

-27 

-25 

a 

0.12 

0.3 

20 

20 

14 

20 

9 

■9 

8 

0.58 

4.6 

113 

-149 

-89 

-50 

9 

0.61 

3.9 

86 

46 

9 

-11 

mSm 

8 

1.69 

0.51 

4.6 

99 

-75 

-50 

-44 

4A 

8 

0.66 

3.7 

119 

24 

-15 

-31 

5A 

8 

1.74 

0.57 

3.7 

91 

-47 

-32 

-12 

6A 

mm 

1.20 

0.60 

5.6 

102 

5 

-38 

-43 

7A 

n 

1.59 

0.50 

4.0 

78 

-59 

-59 

8A 

u 

1.18 

0.50 

4.6 

67 

0 

5 

AVG 

8 

1.45 

0.57 

4.3 

94 

-32 

-34 

-24 

a 

0.05 

0.6 

17 

59 

31 

19 

26 


TABLE  4 


FLUX  SUMMARY  -  CORN  FIELDS  B  AND  C 


Fit 

Run 

Alt 

m 

Dist 

km 

aWge 

mSl 

°co  2 

ppm 

Heat 

Flux 

Watts  m'2 

co2 

Flux  kg  ha'1  h'1 

Trk  2 

Trk  1A 

Trk  IB 

B 

m 

15 

0.66 

3.0 

141 

-27 

-37 

-29 

16 

0.55 

3.9 

-i 

-10 

-9 

17 

0.82 

4.0 

185 

-28 

-31 

-23 

4B 

14 

1.86 

0.50 

4.0 

16 

30 

-8 

1 

AVG 

16 

1.77 

0.63 

3.8 

106 

■a 

-22 

-15 

a 

0.12 

0.4 

64 

mmm 

13 

12 

9 

IB 

9 

5.3 

100 

-3 

-16 

-6 

3B 

7 

1.56 

4.0 

95 

14 

4 

4B 

7 

1.37 

3.5 

22 

-6 

5B 

8 

1.70 

0.62 

3.9 

116 

-5 

iBUr 

6B 

mm 

1.62 

0.49 

4.2 

59 

-15 

mSm 

7B 

mm 

1.49 

0.47 

3.9 

42 

14 

8 

8B 

mm 

1.46 

0.54 

4.2 

53 

-15 

-20 

-11 

AVG 

a 

1.50 

0.59 

4.1 

70 

-2 

-13 

-10 

a 

■ 

0.11 

0.5 

32 

11 

16 

16 

FIELD  C 

m 

4 

1C 

42 

0.63 

3.0 

42 

-57 

-29 

-32 

3C 

21 

0.62 

5.3 

20 

-7 

-16 

-35 

4C 

40 

HU 

0.85 

4.2 

76 

-66 

-65 

-72 

5C 

29 

1.60 

0.62 

3.3 

34 

-26 

-23 

-25 

6C 

28 

1.95 

2.2 

112 

-30 

-46 

-43 

7C 

39 

1.63 

0.82 

3.3 

88 

-23 

-62 

-66 

8C 

20 

1.97 

0.97 

3.5 

210 

-61 

-76 

-98 

9C 

20 

1.56 

0.68 

3.3 

22 

20 

8 

7 

AVG 

30 

1.74 

0.77 

3.5 

76 

-31 

-39 

-46 

a 

0.15 

0.9 

60 

28 

27 

31 

TABLE  5 

FLUX  SUMMARY  -  LAKES 


Fit 

Run 

'  ' 

Alt 

m 

Dist 

km 

ffWge 

mS'1 

aco2 

ppm 

Heat 

Flux 

Watts  m‘2 

C02  FLUX  kg  ha 1  h 1 

Trk  2 

Trk  1A 

Trk  IB 

LAKE  ONTARIO 

5 

4L 

23 

7.53 

0.58 

2.6 

4 

18 

-2 

-2 

5L 

22 

9.99 

0.31 

2.3 

1 

-8 

1 

0 

6L 

33 

6.84 

0.37 

2.1 

1 

1 

0 

3 

7L 

34 

13.33 

0.34 

2.3 

0 

-6 

-8 

-5 

8L 

44 

7.94 

0.42 

2.8 

5 

-3 

-3 

0 

11L 

23 

11.93 

0.36 

2.1 

2 

0 

3 

1 

12L 

38 

8.20 

0.62 

2.8 

-6 

-17 

-5 

-5 

AVG 

31 

9.39 

0.43 

2.4 

1 

to 

-2 

-1 

7 

0.11 

0.3 

3 

10 

3 

3 

BB 

mm 

8 

H 

■ 

2.99 

m 

2.1 

1 

-11 

2 

-3 

0.89 

2.6 

-2 

-7 

1 

1 

K9 

1.93 

la 

2.3 

0 

6 

6 

7 

AVG 

28 

1.93 

0.37 

2.3 

0 

-4 

3 

2 

7 

0.09 

0.2 

1 

7 

2 

4 

28 


TABLE  6 


FLUX  SUMMARY  -  FOREST 


Fit 

Run 

Alt 

m 

Dist 

km 

aWge 

mS1 

ffco2 

ppm 

Heat 

Flux 

Watts  m'2 

C02  Flux  kg  ha'1  h'1 

Trk  2 

TrklA 

Trk  IB 

KINGS' 

K)N  TO  OTTAWA 

6 

IF 

42 

0.86 

4.4 

-6 

-31 

3 

-1 

2F 

36 

0.83 

3.9 

20 

-8 

0 

-5 

3F 

34 

0.87 

3.5 

-9 

24 

7 

-1 

AVG 

37 

10.58 

0.85 

4.0 

2 

-5 

3 

-2 

a 

0.02 

0.4 

13 

23 

3 

2 

LAROSE  FOREST 

7 

IF 

28 

11.88 

0.74 

3.3 

122 

-28 

-37 

-21 

2F 

32 

12.44 

0.90 

3.3 

160 

-23 

-37 

-34 

AVG 

30 

12.16 

0.82 

3.3 

141 

-26 

-37 

-28 

a 

0.08 

0.0 

19 

3 

0 

6 

NEAR] 

HELD  B 

■| 

9 

IF 

23 

0.95 

126 

-47 

-52 

-42 

2F 

28 

0.85 

94 

50 

18 

11 

3F 

25 

1 

0.83 

135 

-65 

-71 

-58 

4F 

18 

0.71 

0.65 

53 

14 

14 

7 

5F 

19 

1.56 

fpJIiUMt 

260 

-5 

-43 

-45 

7F 

18 

1.82 

0.93 

4.0 

157 

-19 

0 

-2 

8F 

23 

1.37 

0.92 

4.0 

106 

34 

67 

68 

AVG 

22 

1.38 

0.88 

3.7 

134 

-5 

-10 

-9 

a 

0.11 

0.5 

60 

39 

45 

40 

TABLE  7 


FLUX  SUMMARY  -  RICHMOND  FARMLAND,  MER  BLEUE 


Fit 

Run 

Alt 

m 

Dist 

km 

0Wge 

m  S'1 

aco2 

ppm 

Heat 

Flux 

Watts  m'2 

co2 

Flux  kg  ha'1  h'1 

Trk  2 

Trk  1A 

Trk  IB 

RICHMOND  FARMLA 

ND 

8 

1R 

29 

5.55 

0.64 

3.0 

61 

-28 

-25 

2R 

26 

5.65 

0.82 

3.3 

85 

-70 

-32 

3R 

29 

3.29 

0.76 

3.2 

82 

-21 

-18 

4R 

26 

3.06 

0.62 

3.5 

61 

-9 

-23 

-5 

5R 

32 

6.29 

0.62 

3.2 

78 

-36 

-40 

-17 

AVG 

28 

4.77 

0.69 

3.2 

73 

-33 

-30 

-19 

a 

0.08 

0.2 

10 

21 

9 

9 

MER  BLEUE 

2M 

18 

5.94 

0.61 

3.0 

95 

-10 

-12 

-11 

(1)  At  flight  level 

(2)  Data  in  last  4  columns  from  Tables  3  to  7,  computed  from  spectra  and  cospectra  10  -  800  m  band. 
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FIG.  2:  AGRICULTURE  CANADA  OPEN-PATH  CO.  ANALYZER 


CF-POK  AGRI  C02  1982 
RUN  04-08C  CORN  FIELD  C 
RECORD  DATE  AUG23/82 

ANALYSIS  DATE  0CT25/82 
PLOT  DATE  22/02/83 


TIME  0.0 


5 


10 


15 


2 


3 


IME  HISTORIES,  LAROSE  FOREST,  RUN  07-011 


CF-POK  AGRI  C02  1982 
RUN  07-02M  MER  BLEU 
RECORD  DATE  AUG31/82 

ANALYSIS  DATE  0CT19/82 
PLOT  DATE  22/02/83 


TIME  0.0 


AVERAGE  COSPECTRA  -  CORN  FIELD  B 


KKS 


FIG.  12:  AVERAGE  COSPECTRA  -  LAKES 


FIG.  13:  AVERAGE  COSPECTRA  -  LAROSE  FOREST 


LT.  22  i 
2.9 


AVERAGE  COSPECTRA  -  FOREST  NEAR  B 


FIG.  15:  AVERAGE  COSPECTRA  -  ALL  FOREST  RUNS 


FIG.  16:  AVERAGE  COSPECTRA  •  RICHMOND  FARMLAND 


rjt 


FIG.  19:  TIME  HISTORIES  FOR  THREE  RUNS  USED  FOR  SIDESLIP  STUDIES 
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SUMMARY/SO  MM  AIR  E 

At  the  request  of  Agriculture  Canada,  the  Flight  Research 
Laboratory  of  the  National  Aeronautical  Establishment  (NAE)  has  operated 
the  Twin  Otter  atmospheric  research  aircraft  fitted  with  a  C02  analyzer 
in  order  to  assess  the  feasibility  of  making  airborne  carbon  dioxide  flux 
measurements.  This  report  will  summarize  progress  to  date  and  present  an 
analysis  of  the  data  collected  during  a  series  of  nine  test  flights  in  August 
and  September,  1982. 

The  experimental  results  show  that  it  is  possible  to  make  realistic 
measurements  of  the  heat  and  C02  fluxes  from  an  aircraft  instrumented  to 
measure  vertical  gusts.  Special  tests  and  analysis  procedures  were  conducted 
to  investigate  the  accuracy  of  the  carbon  dioxide  measurements.  Recom¬ 
mendations  for  improvements  in  the  design  of  future  teBt  programs  are 
made, 
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